Large-scale application of free energy perturbation calculations for antibody design.

Alchemical free energy perturbation (FEP) is a rigorous and powerful technique to calculate the free energy difference between distinct chemical systems. Here we report our implementation of automated large-scale FEP calculations, using the Amber software package, to facilitate antibody design and evaluation. In combination with Hamiltonian replica exchange, our FEP simulations aim to predict the effect of mutations on both the binding affinity and the structural stability. Importantly, we incorporate multiple strategies to faithfully estimate the statistical uncertainties in the FEP results. As a case study, we apply our protocols to systematically evaluate variants of the m396 antibody for their conformational stability and their binding affinity to the spike proteins of SARS-CoV-1 and SARS-CoV-2. By properly adjusting relevant parameters, the particle collapse problems in the FEP simulations are avoided. Furthermore, large statistical errors in a small fraction of the FEP calculations are effectively reduced by extending the sampling, such that acceptable statistical uncertainties are achieved for the vast majority of the cases with a modest total computational cost. Finally, our predicted conformational stability for the m396 variants is qualitatively consistent with the experimentally measured melting temperatures. Our work thus demonstrates the applicability of FEP in computational antibody design.

Development of potent and effective synthetic SARS-CoV-2 neutralizing nanobodies.

The respiratory virus responsible for coronavirus disease 2019 (COVID-19), severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has affected nearly every aspect of life worldwide, claiming the lives of over 3.9 million people globally, at the time of this publication. Neutralizing humanized nanobody (VHH)-based antibodies (VHH-huFc) represent a promising therapeutic intervention strategy to address the current SARS-CoV-2 pandemic and provide a powerful toolkit to address future virus outbreaks. Using a synthetic, high-diversity VHH bacteriophage library, several potent neutralizing VHH-huFc antibodies were identified and evaluated for their capacity to tightly bind to the SARS-CoV-2 receptor-binding domain, to prevent binding of SARS-CoV-2 spike (S) to the cellular receptor angiotensin-converting enzyme 2, and to neutralize viral infection. Preliminary preclinical evaluation of multiple VHH-huFc antibody candidates demonstrate that they are prophylactically and therapeutically effective in vivo against wildtype SARS-CoV-2. The identified and characterized VHH-huFc antibodies described herein represent viable candidates for further preclinical evaluation and another tool to add to our therapeutic arsenal to address the COVID-19 pandemic.

Discovery of Small-Molecule Inhibitors of SARS-CoV-2 Proteins Using a Computational and Experimental Pipeline.

A rapid response is necessary to contain emergent biological outbreaks before they can become pandemics. The novel coronavirus (SARS-CoV-2) that causes COVID-19 was first reported in December of 2019 in Wuhan, China and reached most corners of the globe in less than two months. In just over a year since the initial infections, COVID-19 infected almost 100 million people worldwide. Although similar to SARS-CoV and MERS-CoV, SARS-CoV-2 has resisted treatments that are effective against other coronaviruses. Crystal structures of two SARS-CoV-2 proteins, spike protein and main protease, have been reported and can serve as targets for studies in neutralizing this threat. We have employed molecular docking, molecular dynamics simulations, and machine learning to identify from a library of 26 million molecules possible candidate compounds that may attenuate or neutralize the effects of this virus. The viability of selected candidate compounds against SARS-CoV-2 was determined experimentally by biolayer interferometry and FRET-based activity protein assays along with virus-based assays. In the pseudovirus assay, imatinib and lapatinib had IC50 values below 10 µM, while candesartan cilexetil had an IC50 value of approximately 67 µM against Mpro in a FRET-based activity assay. Comparatively, candesartan cilexetil had the highest selectivity index of all compounds tested as its half-maximal cytotoxicity concentration 50 (CC50) value was the only one greater than the limit of the assay (>100 µM).